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Abstract 

The first result of the pp2pp experiment at RHIC on elastic scattering of polarized protons 
at \fs = 200 GeV is reported here. The exponential slope parameter b of the difFractive peak 
of the elastic cross section in the t range 0.010 < \t\ < 0.019 (GeV/c) 2 was measured to be 
b = 16.3 ±1.6 (stat.) ± 0.9 (syst.) (GeV/c)- 2 . 



PACS numbers: 13.75Cs, 29.27Hj, 14.20Dh 



Although elastic scattering has been measured in pp collisions up to y/s = 1.8 TeV, the 
highest energy pp data reach only to 63 GeV. We present here the first measurement of 
the slope parameter b in forward proton-proton elastic scattering obtained by the pp2pp 
experiment at the Relativistic Heavy Ion Collider (RHIC) at ^/s = 200 GeV. 

The pp2pp experiment |lj is designed to measure polarized pp elastic scattering at 
RHIC, which will provide proton beams with polarizations of 0.7 and luminosities up to 
2 x 10 32 cm~ 2 sec _1 . The main goal of the experiment is to study the spin dependence of 
elastic scattering in the squared four-momentum transfer range 4x 10 -4 < \t\ < 1.3 (GeV/c) 2 
and 50 < y/s < 500 GeV. 

By measuring elastic scattering of polarized protons in the nonperturbative regime of 
QCD at RHIC, one has a unique opportunity to probe the spin structure of the nucleon 
and of the exchanged mediators of the force, the Pomeron and its odd C-parity partner, the 
Odderon. The pp2pp experiment, part of the RHIC spin program, studies the physics of 
elastic scattering and diffractive dissociation. It addresses the main unsolved problems in 
particle physics- long range QCD and confinement. 

The slope b for \t\ < 0.5 (GeV/c) 2 is inherently sensitive to the exchange process, and its 
dependence on y/s will allow to distinguish among various QCD based models of hadronic 
interactions. Some interesting features of b observed in pp are not yet confirmed in pp elastic 
scattering. In general, the forward peak does not show a simple exponential behavior. The 
t distribution becomes less steep as \t\ increases from 0.02 (GeV/c) 2 to 0.20 (GeV/c) 2 , 
although at the highest Tevatron energies this was not observed. It is therefore of interest 
to see the b behavior in the RHIC energy range. 

In RHIC the two protons collide at the interaction point (IP), and since the scattering 
angles are small, scattered protons stay within the beam pipe of the accelerator. They follow 
trajectories determined by the accelerator magnets until they reach the detectors, which 
measure the x, y coordinates in the plane perpendicular to the beam axis. The coordinates 
are related by the beam transport equations to the corresponding quantities at the IP: 

x = an ■ x + L x eff ■ 9* x + a 13 ■ y + a u ■ 9* 

y = a 31 -x + a 32 ■ 9* x + a 33 • y + L v eff ■ 9* (1) 

where xo, yo and 9*, 9* are the positions and scattering angles at the IP and a»j and L e // 
are the elements of the transport matrix. The optimum condition for the experiment is 
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FIG. 1: Layout of the pp2pp experiment. Note the detector pairs RP and RP' lie in different 
RHIC rings. Scattering is detected in either one of two arms: Arm A is formed from the upper 
half of RP' and the lower half of RP. Conversely, Arm B is formed from the lower half of RP' and 
the upper half of RP. 

to minimize the dependence of the measured coordinates on the unknown collision vertex, 
i. e. to have the o^-'s small and the L e //'s as large as possible. In that case, called "parallel 
to point focusing", rays that are parallel to each other at the interaction point are focused 
nearly to a single point at the detector. Since in practice such a condition is achieved for 
one coordinate only, in our case y, Eq. (JTJ) then simplifies to y « L v e ^ • 8*. 

The momentum-transfer interval for the data presented here is 0.004 < |£| < 0.032 
(GeV/c) 2 . In our 14 hour run of January 2002, the RHIC orbit betatron function [J at the 
IP was (3* = 10 m, resulting in L y e ^ « 24 m. At larger momentum transfers the acceptance 
is limited by the aperture of the RHIC focusing quadrupoles. 

The layout of the experiment is shown in Fig. 1. The identification of elastic events 
is based on the collinearity criterion, hence it requires the simultaneous detection of the 
scattered protons in the pair of Roman Pot (RP) detectors [3], RP and RP', on either side 
of the IP. Additionally, a set of scintillators located outside of the beam pipe near the IP 
provide detection of inelastic events. 

The RP's are insertion devices allowing four silicon strip detectors (SSD) to be positioned 
just above and below the beam orbits. The SSD's inside the pots record the x, y coordinates 



of the scattered protons. The silicon detectors are made of 0.40 mm thick n-type silicon 
with p + -type implanted strips of 0.07 mm width and a strip pitch of 0.10 mm. Two of the 
detectors have 512 strips implanted along the longer side of the rectangle, the other two 768 
strips perpendicular, resulting in an active area of 75 x 45 mm 2 . Each strip is capacitively 
coupled to an input channel of a SVXIIe jj], which has 128 channels with preamplification, 
a 32 event pipeline, and a Wilkinson-type ADC. 

The amount of charge collected due to a 100 GeV/c proton passing through the silicon 
detector corresponds to an energy deposit of about 200 keV. In 80% of the events, this 
deposited energy is confined to a single strip, and otherwise shared between neighboring 
strips if the particle passed through a 30/im wide region in between the strips. 

The elastic trigger scintillators were 8 mm thick, 80 x 50 mm 2 in area, and were viewed by 
two photomultiplier tubes. To produce a highly efficient and uniform trigger the two signals 
from the tubes formed a logical OR. The elastic event trigger is a coincidence between 
signals in the RP's scintillators, belonging either to arm A or arm B (see Fig. 1). The 
trigger efficiency was greater than 0.99. For each event, time and amplitude were digitized 
and recorded. 

The coordinate in the SSD is calculated as an energy-weighted average of the positions 
of the hit strips. Clusters of more than three hit strips were excluded. The detection 
efficiency for every SSD strip was calculated using the redundancy of the silicon planes for 
identification of elastic events. The average silicon detector plane efficiency for arm A was 
0.97. 

The collinearity of elastic events implies that the two coordinates obtained from the 
silicon detectors on either side of the interaction point are correlated. This correlation is 
shown for the y coordinates from arm A in Fig. |21 The widths of the coordinate difference 
distributions, a x and o y , were determined. Events for which a/ Ax 2 + Ay 2 < 4 \/o 2 x + o" 2 
were retained for the analysis. The widths are dominated by the beam angular emittance 
of about 127T /im and by the uncertainty of about 60 cm (rms) in the vertex position along 
the beam axis. 

At least six of the possible eight planes were required to have hits to be accepted for 
elastic events. Out of 196,000 elastic triggers for arm A about 84% were reconstructed. Most 
reconstruction failures are accounted for by the larger area of the scintillator compared to 
the active area of the SSD packages. The above mentioned correlation cut of 4cr removed 
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FIG. 2: Correlation between the y coordinates as measured by the two detectors of arm A for 
elastic events before cuts being applied. 

another 3.8%, while the requirement of six hit planes contributing to the track reconstruction 
cut another 0.3%. To reduce the contamination of the elastic event sample with tracks from 
background particles, not more than two planes with more than one hit per event were 
accepted. This reduced the event sample by another 3.2%, giving a total of 153,000 elastic 
events for this arm. A similar analysis was carried out for arm B, but because of the noise 
level being considerably higher, it was used only for consistency checks, but not included in 
the final analysis presented here. 

For each event the scattering angle 9 and azimuth <ft were calculated for each proton and 
then averaged. The scattering angle is related to the square of the four-momentum transfer, 
t, via —t^(p-8) 2 . A restriction of the <fi range leads to a uniform geometric acceptance in 
a limited t-range. For 45° < <p < 135° or 225° < <p < 315° that range is 0.010 < |t| <0.019 
(GeV/c) 2 . In Fig. |3]the correlation between t and <fi is shown for reconstructed events. The 
determination of the slope parameter b is confined to the t region for which no acceptance 
correction is required. The final selection therefore yields 58,511 events. The uncorrected 
dN/dt distribution resulting from the 0-cut is shown in Fig. Altogether with the acceptance 
function obtained from Monte Carlo studies. 
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FIG. 3: Correlation between t and <b for reconstructed events. 



The differential cross section d<r/dt for elastic scattering in the forward angle region is 
determined by Coulomb and nuclear amplitudes and the interference term between them. 
The cross section is given by (see for example Ref. Jl]]) 



da „ faG 

— = An (he) 
dt K ' 



,;2 



+ 1+P o 2 -e-^l 
16 tt {hef tot 



- ( P + A^).^.a tat -e-^\ (2) 
with a the fine structure constant, Ge the electric form factor of the proton, A$ the Coulomb 



phase 6J, p the ratio of the real to imaginary part of the forward scattering amplitude, o tot 
the total cross section, and h the nuclear slope parameter. The dominant contribution in 
our t region is the second term in this expression. 

A least squares fit was performed to the distribution of Fig. 0] using Eq. (J2J) with h and 
a normalization constant as free parameters. The total cross section and p were fixed to 
o t ot — 51-6 mb [3] and p = 0.13 8]. These values of a tot and p come from fits to the existing 
pp data taken at energies below 63 GeV and world pp data. 
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FIG. 4: The distribution of dN/dt within the (ft region selected as described in the text. The 
two distributions shown are the measured data and the simulated acceptance function. The fit is 
shown by the solid line. 



The resulting slope parameter is 

b = 16.3 ± 1.6 (stat.) ± 0.9 (syst.) (GeV/c)" 2 . 

The evaluation of the systematic errors due to the uncertainty in beam emittance, vertex 
positions and spread, beam transport matrix elements, and incoming beam angles was based 
on Monte Carlo simulations. These simulations used the geometry of the experimental setup 
and efficiency of the detectors as an input. The largest single source of the systematic error 
was the uncertainty of the initial colliding beam angles. 

There is also a correlation between b and the values of p and a to t- We find that the 
changes Ap = ±0.02 and Act = ±4 mb result in changes in b of Ab = T-0.32 (GeV/c) -2 and 
Ab = T-0.07 (GeV/c)- 2 , respectively. 

An independent analysis of the data was performed using different selections of hits and 
elastic events. In particular, a t-dependent cut on eft was applied, which allowed an increase 
in the t range and the number of accepted elastic events. The b slope values obtained from 
both analyses agree within statistical errors. 

Our result for the slope parameter b is shown in Fig. El together with the world data on 




FIG. 5: The result for the slope parameter b of this experiment compared to the world pp and pp 
data set. The data are drawn from the Durham Database Group (UK). Only statistical errors are 
shown. 

elastic pp and pp scattering. This result is about one standard deviation higher than an 
extrapolation of world data to the energy of this experiment Q], 0, [lfll ] . 

In the future, a full complement of two sets of Roman Pot detector pairs will be used, 
two pairs at each side of the IP, to allow a direct measurement of the scattering angles. This 
will reduce the systematic error due to the uncertainty of the vertex position. An expected 
increase of the RHIC luminosity will result in a reduction of the statistical error and will 
make the studies of the polarized observables An and Ann feasible. 
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